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We study the contribution to the entanglement entropy of (2+l)-dimensional conformal field 
theories coming from a sharp corner in the entangling surface. This contribution is encoded in a 
function a{9) of the corner opening angle, and was recently proposed as a measure of the degrees of 
freedom in the underlying CFT. We show that the ratio a{9)/CT, where Ct is the central charge in 
the stress tensor correlator, is an almost universal quantity for a broad class of theories including 
various higher-curvature holographic models, free scalars and fermions, and Wilson-Fisher fixed 
points of the 0{N) models with N = 1,2,3. Strikingly, the agreement between these different 
theories becomes exact in the limit 9^ tv, where the entangling surface approaches a smooth curve. 
We thus conjecture that the corresponding ratio is universal for general CFTs in three dimensions. 


Many interacting gapless quantum systems do not pos¬ 
sess simple particle-like excitations, making it difficult 
to quantify their effective number of degrees of freedom 
(dof) at low-energy. Conformal field theories (CFTs) 
constitute an important example. For CFTs in two 
spacetime dimensions (2d), the Virasoro central charge 
is a good measure of the dof. It appears in many quan¬ 
tities, such as the thermal free energy and the entangle¬ 
ment entropy (EE), and decreases under renormalization 
group (RG) flow [1]. In higher dimensions, the concept of 
quantum entanglement is emerging as a fundamental di¬ 
agnostic for such measures [2, 3]. E.g., it was instrumen¬ 
tal in finding an analogous RG monotone for 3d CFTs, 
with the EE of a disk-shaped region [4]. We shall study 
another measure of recent interest [5-16]: the coefficient 
capturing the contribution of sharp corners to spatial en¬ 
tanglement. 

In the context of quantum field theory, the EE is de¬ 
fined for a spatial region E as: S = —Tr (py In py), where 
pv is the reduced density matrix produced by integrat¬ 
ing out the dof in the complementary region V. In the 
groundstate of a 3d CFT, the EE takes the form: 

S = B£/S-a{e)\n{e/6) + 0{l), (1) 

where 5 is a short-distance cutoff, e.g., the lattice spacing, 
and £, a length scale associated with the size of V. The 
first, ‘area law’, term depends on the UV regulator and 
scales with the size of the boundary. The second one 
appears only when V has a sharp corner with opening 
angle 0 G [0, 27r), Fig. 1. Crucially, a{9) is a regulator 
independent coefficient that characterizes the underlying 
CFT. It is positive and satisfies a{2'K — 9) = a{9) [5], and 
behaves as follows: 

a(0 —>■ tt) ~ (T (tt — 0)^ ; a(0 —>■ 0) ~ k/0 (2) 

in the limits of a nearly smooth entangling surface and 
a very sharp corner, respectively. It has been studied 
for a variety of systems: free scalars and fermions [5-7], 
interacting scalar theories via numerical simulations [8- 
10], Lifshitz quantum critical points [11], and holographic 
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FIG. 1: a) An entangling region V of size £ with a corner; 
b) The holographic entangling surface 7 for a region on the 
boundary of AdS 4 with a corner. 




models [12]. The results suggest that a{9) is an effective 
measure of the dof in the underlying CFT [7, 10]. 

Another quantity measuring dof is the central charge 
Ct, associated with the stress tensor of the CFT. It 
characterizes the vacuum two-point function: 

{T)_i,^{x)Txp{0)) = ir^j/,Ap(a;), (3) 

where Ip,v,\p is a dimensionless tensor structure fixed by 
symmetry [17]. In the following, we will show that the 
ratio a{9)/CT is almost universal for a broad class of 
theories. In fact, we find that this agreement becomes 
exact in the limit 0 —>■ tt. Hence, using (2) we conjecture 
there is a universal ratio, i.e., the same value arises in 
any general 3d CFT: 

a/CT = 7rV24 . (4) 


This is a striking result since the EE can generally be 
regarded as a nonlocal quantity but our analysis indicates 
that the regulator independent corner contribution to the 
EE is controlled by a local correlation function (3). 

Holographic calculations: The AdS/CFT cor¬ 
respondence posits that the physics of certain d- 
dimensional CFTs has an equivalent description in terms 
of gravity coupled to a negative cosmological constant in 
d-£l dimensions. In such holographic CFTs, EE is com¬ 
puted using the Ryu-Takayanagi prescription [18, 19]: 


S{V) = ext 

‘■yr-.jV 


AG 


( 5 ) 
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That is, given a region V in the boundary CFT, we 
consider all codimension -2 surfaces 7 in the dual AdS 
spacetime which are homologous to V on the asymp¬ 
totic boundary and find the surface which extremizes the 
area functional within this class. The EE is then given 
by Bekenstein-Hawking formula A( 7 )/( 4 G') evaluated on 
this extremal surface, where G is the gravitational con¬ 
stant, Fig. lb. This prescription for the holographic EE 
(HEE) has been used to compute the corner coefficient 
aTs{9) for 3d CFTs dual to 4d AdS [12, 20]; where ‘E’ 
indicates that the bulk theory is described by Einstein 
gravity. a^{ 6 ) is only implicitly known in terms of 2 inte¬ 
grals (appendix A) but it is easily evaluated numerically 
and the result is plotted in Fig. 2. Evaluating Oe^O) for 
the limits in ( 2 ), we find: = r(|)‘^ Z^/(27rG) and 

(Te = T^/( 87 rG), where L is the radius of curvature of 
the dual AdS geometry [ 21 ]. 

Just as for and aE(0) carries an overall factor of 
the ratio L^/G which is indicative of the number of dof in 
the boundary CFT. E.g., [22] provides a top-down model 
where the boundary CFT is a specific supersymmetric 
gauge theory and L^/G ~ A^/-\/A where N and A are 
the rank of the gauge group and the ’t Hooft coupling, 
respectively. However, this same ratio appears in many 
other physical quantities, e.g., the stress tensor two-point 
function (3), EEs for general regions or the thermal en¬ 
tropy density. Essentially, this ubiquitous appearance of 
I? jG occurs because it is the only dimensionless param¬ 
eter in the bulk Einstein theory. As a result, all of these 
a priori distinct measures of the dof are all related in this 
class of CFTs. However, by introducing higher curvature 
interactions, the bulk gravity acquires new dimensionless 
couplings and each of the above measures can acquire a 
distinct dependence on these couplings. Hence we can 
see whether the various measures are independent or if 
they encode the same information, e.g., see [23-25]. 

Higher curvature interactions appear generically in 
string theoretic models, e.g., as o' corrections in the low- 
energy effective action [26]. However, rather than con¬ 
structing explicit top-down holographic models, our ap¬ 
proach will be to examine simple toy holographic models 
involving higher curvature gravity in the bulk. Our per¬ 
spective is that if there are interesting universal proper¬ 
ties which hold for all CFTs, then they should also appear 
in the holographic CFTs defined by these toy models as 
well. This approach has been successfully applied before, 
e.g., in the discovery of the F-theorem [2, 3]. 

In particular, we focus our attention on the following 
simple gravitational theory 


1 = 


\ 6 

IGttG L2 


+ R + L^XiRX^, + 


( 6 ) 


The first two terms above are the cosmological constant 
with A = and the standard Einstein term. The 

next two interactions are controlled by the dimension¬ 
less couplings Ai ,2 and contain X 4 = R^upaR^^^^ — 
AR^yR^^ + R^, which is the Euler density on 4d man¬ 
ifolds. While X 4 alone would be topological, RX 4 and 


X 4 are not, and so these terms do modify the gravi¬ 
tational equations of motion. In particular, AdS space 
(with radius of curvature L) is a solution provided: 
1 = L'^/P - 24AiLVF® + 96 A 2 TVT®. Note that in 
the following, we treat the higher curvature interactions 
as perturbative corrections and we will only calculate to 
leading order in A 1 . 2 . 

With higher curvature interactions, the Ryu- 
Takayanagi prescription (5) must be modified [27-31]. 
Schematically, we have S{V) = extS'g„v( 7 ) where the 

new entropy functional contains the ‘higher curvature 
corrections’ to the Bekenstein-Hawking formula. For the 
above action (A 8 ), one finds [ 21 , 32] 




rdPyVh 
. 4G 


l-fAlL'‘(A’4-f2i?7^)-f4A2T®A’47^ , (7) 


where h^j and TZ are the induced metric and intrinsic 
Ricci scalar, respectively, on 7 . 

When evaluating the HEE in the ground state of the 
boundary CFT, the bulk geometry is empty AdS 4 , and 
(A9) simplifies because X 4 = 2A/L‘^ and R = —Vljl? are 
constants. Further TZ is topological on two-dimensional 
manifolds, i.e., 'JhTZ can be written as a total derivative. 
Thus, in this case, the HEE is still determined by an ex¬ 
tremal area surface in the bulk and the corner term con¬ 
tribution is easily evaluated for the CFTs dual to (A 8 ): 

a{9) = aaE{d) with a = 1-|-24Ai-I-G(A^). ( 8 ) 


Hence the corner coefficient is modified by an overall fac¬ 
tor but the 0-dependence is unchanged. The above result 
also implies k, = a Ke and cr = a ctb • 

Now we compare this result to other measures of the 
dof in the holographic CFTs dual to the action (A 8 ), 
beginning with the central charge G^,. Since the stress 
tensor in the boundary CFT is dual to the metric per¬ 
turbations in the bulk gravity theory, (3) translates to 
a statement about the graviton propagator between two 
boundary points in AdS 4 [24, 33]. Hence we must deter¬ 
mine the normalization of the graviton kinetic term in 
our higher curvature theory (A 8 ). However, we should 
note that generically the higher curvature interactions 
introduce additional massive dof in the bulk theory with 
~ l/(AiL^) [3, 21]. Hence, the metric will couple to 
the stress tensor, but also to some additional (nonuni¬ 
tary) operators. However, in a perturbative framework, 
it is straightforward to determine Gt by computing the 
equation of motion corresponding to the massless spin -2 
graviton mode. In particular, we consider metric fluc¬ 
tuations g^i, = g^y h^y, where is the background 
AdS 4 metric and h^y 1. The equation of the physical 
massless graviton can be isolated by choosing the trans¬ 
verse gauge condition W^h^y = yh, supplemented by 
the tracelessness condition g^'^h^y = 0 , and reads 

h^y=S'KGT^y, (9) 
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where a is the same coefficient as in ( 8 ). We have in¬ 
cluded the bulk stress tensor of matter fields to es¬ 
tablish the normalization of the graviton kinetic term. 
The net effect of the appearance of a in (9) is to modify 
the holographic result for the stress tensor correlator (3) 
by an overall factor a: 

C, = (l + 24Ai+0(A?))C'^.,, (10) 

where the Einstein result is = 3 L^/( 7 r^G) [24, 33]. 

Remarkably then, for all the holographic CFTs dual 
to (A 8 ), we find that a{9)/CT = aE{0)/CT,E, i-C., we find 
a universal ratio that is independent of the details of 
the theory. It is notable that this universality does not 
occur when considering other measures of the boundary 
dof. Here we examined the coefficient in the thermal 
entropy density, s = CgT^, which can be determined by 
evaluating the horizon entropy of a planar AdS black 
hole, and the RG monotone E, which is determined using 
the HEE prescription with a circular entangling surface 
in the boundary. The results are 

Cs = (l + 24Ai -|- I 6 A 2 -f 0(A^)) Cs_E , 

F = (l-h48Ai-96A2-bG(A2)) Fb, (11) 

where Cs,E = 47 r^L^/( 9 G) and Ee = E^/(2G). Hence each 
of these measures of dof has a unique signature in terms 
of the couplings Ai ,2 and it is only the ratio a{0)/CT that 
yields a universal result. Further in [21], this holographic 
analysis was extended to an 8 -parameter family of higher 
curvature theories and the same universal ratio arises in 
the dual boundary CFTs. 

QFT comparison: The universality revealed by our 
holographic analysis suggests more broadly that a{9)/CT 
provides a useful normalization if we wish to compare the 
corner coefficients of different field theories. The only 
cases where a{9) is known for a wide range of angles are 
the free massless scalar and fermion [5-7]. Despite be¬ 
ing free field theories, the required calculations are tech¬ 
nically demanding and a{9) is only given implicitly in 
terms of a complicated set of nonlinear differential and 
algebraic equations. Explicit results were given for Tay¬ 
lor expansions of a{9) around = tt (to 14th order), 
and for 9 = 7 r/ 4 , 7 r/ 2 , 37 r /4 [7], which were obtained us¬ 
ing the numerical methods of [34]. The values of k, i.e., 
the coefficient at 0 —)■ 0 , can also be determined using 
a conformal transformation [ 6 , 21]. Finally analytic ex¬ 
pressions for the free field central charges are known [17]: 
(^scalar ^ 3/(327r2) and = 3 /( 167 r 2 ). Combin¬ 

ing these results, the ratios a{9)/CT for the holographic 
and free field theories are compared in Fig. 2. The solid 
lines for the free fields are interpolating functions which 
combine the Taylor expansions around 9 = tt with the 
coefficients k a.t 9 = 0. 

Fig. 2 also includes a{'K 12 )ICt for the Wilson-Fisher 
fixed points of the 0{N) models with = 1,2, 3. In this 
case, a(7r/2) was evaluated numerically using state of the 
art numerical simulations for lattice Hamiltonians with 
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FIG. 2: (Top) aid) ICt from holography (gray), a free Dirac 
fermion (red) and a scalar (blue), plus the corresponding lat¬ 
tice data points obtained numerically [7] (red/blue squares). 
We also show a{'K 12) jC t for the E = 1,2,3 Wilson-Fisher 
0{N) CFTs, and the trial function (14) (purple). (Bottom) 
Same quantities normalized by \a{6 )/Ct]]ioIo- 

the corresponding quantum critical points [8-10]. Next, 
Ct was recently determined with great accuracy using 
conformal bootstrap methods [35]. Remarkably the lat¬ 
tice results indicated that a(7r/2) satisfies the emergent 
scaling {n/ 2 ) ~ iVa^®‘"®(7r/2), to within the numer¬ 

ical accuracy. Further, the bootstrap calculations yield 
Ct with the same approximate scaling for 1 < < 3. 

Hence a{'K 12 )JCt is nearly independent of N for these 
CFTs, Table I. Finally, Fig. 2 also shows a trial function, 
with a simple closed form, obtained from [36-38], see 
(14). 

We see that all the results in Fig. 2 are in remarkable 
agreement for the whole range of 9. Refs. [7, 39] had 
previously noted a qualitative similarity between the cor¬ 
ner coefficient in the holographic and free field theories. 
However, here we see a remarkable quantitative agree¬ 
ment after normalizing a{9) with the central charge Ct- 
The free scalar curve differs from the holographic func¬ 
tion by no more than 13% in the whole range, while the 
agreement with the free fermion is even better, differing 
from the holographic result by no more than 2.5%. This 
excellent agreement also extends to the 0{N) Wilson- 
Fisher CFTs. As we can see in Table I and the lower 
plot in Fig. 2, the ratio a(7r/2)/GT for these CFTs agrees 
remarkably well amongst each other and deviates only 
slightly from the holographic and free field results. 

A conspicuous feature shown in Fig. 2 is that the ra¬ 
tio [a{9 )/Ct\ / [a{9)/CT\hoio is a monotonically decreas¬ 
ing function of the opening angle for both free field the¬ 
ories. Hence the maximum discrepancy between these 
theories occurs ai 9 — Q. The most striking aspect of the 
plot is that these ratios approach 1 as 0 —tt, indicating 
the possibility of universal behavior that extends beyond 
holography to general 3d CFTs. 
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Ising {N=l) 

XY {N = 2) 

Heisenberg (A = 3) 

Free scalar 

Dirac fermion 

AdS/CFT 

a{'Kl2) ICt 
a2(7r/2)/CT 

1.36(14) [S. 35] 
0.62(6) [8, 13, 14, 35, 40, 41] 

1.3(1) [ 9 , 35 ] 
0.62(6) [ 9 , 15, 35 , 42] 

1.3(1) [ 10 , 35 ] 
0.61(6) [ 10 , 16 , 35 , 42] 

1.245 [7, 17] 
0.674 [7, 17] 

1.226 [6, 17 ] 

1.222 [ 12 , 33 ] 


TABLE I: Ratio an{'K 12) jCr for n = 1,2 in different critical theories. The first 3 are the 0{N) Wilson-Fisher CFTs. 


Universal ratio: As described in (2), a{9) is con¬ 
strained to vanish quadratically at d = tt, i.e., in the 
limit of a smooth entangling surface. Our holographic 
analysis indicates that the ratio of the corresponding co¬ 
efficient with the central charge is universal for a broad 
family of holographic CFTs (see also [21]), 

a/C^ = ct^/Ct,e = 7rV24 ~ 0.411234. (12) 

Fig. 2 suggests that the free scalar and free fermion the¬ 
ories yield the same ratio and so we can evaluate this 
ratio for these theories, using the values of Cj- given 
above and those of a given in [5-7]: cTscaiar — 0.0039063, 
CTfermion ^ 0.0078125. We find: (a/Cr)scalar 0.411235, 
(ct/Ct) fermion — 0.411234. Hence the free field ratios 
agree with the holographic result to at least five signif¬ 
icant figures. Recall that the free field values of Ct are 
exact but the values of a are approximate numerical re¬ 
sults [5-7] . Therefore the precision of the agreement here 
is as good as possible. 

This motivated the conjecture (4) that the ratio 
ct/Ct = 7r^/24 for all 3d CFTs. Our conjecture can be 
used to predict the exact free field values for a: 

i^scalar — 1/256 , CTfermion — 1/128 . (1^) 

To test this prediction, we revisited the original compu¬ 
tations for these coefficients [5-7]. In appendix B, we ex¬ 
plicitly give the monstrous integrals needed to calculate 
a for the free fields. Improving the numerical accuracy 
in evaluating these integrals, we found that the agree¬ 
ment between the numerical results and (13) was easily 
extended to 1 part in 10^^ - the accuracy to which we 
limited ourselves. [57] 

These two integrals are very complicated and not sim¬ 
ilar. The fact that they yield simple fractions and differ 
by a factor 2 is highly non-trivial. 

Conclusions & outlook: We have shown that 
a{9)/CT is an almost universal function of the open¬ 
ing angle for a broad class of CFTs, including a fam¬ 
ily of higher-curvature holographic models, free massless 
scalars and fermions, and Wilson-Fisher fixed points of 
the 0{N) models with N = 1,2,3 - see Fig. 2. A striking 
aspect of this result is that generally the EE is consid¬ 
ered a nonlocal quantity but here we found the regulator 
independent corner contribution to the EE encodes es¬ 
sentially the same counting of dof as the local correlation 
function (3). [58] 

While the ratios a{9)/CT for the different theories 
agree well across the whole range of 9, the agreement 
becomes exact as 6* —>■ tt. Given the extremely dif¬ 
ferent nature of these CFTs and the computations in¬ 
volved in evaluating a{9) in each case, we conjectured 


that a/Cr = 7r^/24 is a universal result for all 3d CFTs. 
Checking its validity in additional theories is an exciting 
issue to address. It seems possible to prove our conjec¬ 
ture [43] using the techniques developed in [44, 45]. Fur¬ 
ther, holographic calculations suggest that similar uni¬ 
versal behavior may also arise in higher dimensions [43]. 

Higher Renyi entropies also contain regulator indepen¬ 
dent corner functions a„(0), analogous to a{9) (which 
corresponds to n = 1). With results from lattice stud¬ 
ies of the Wilson-Fisher CFTs [8-10] and calculations 
for the free scalar field [7], Table I suggests that the al¬ 
most universal behavior found for a{9)/CT may extend to 
0'n{9) / Ct with n ^ 1. Further, at d = tt, the Renyi corner 
terms also have a quadratic zero, an{9 ^ tt) = (j„(7r —0)^. 
It would be interesting to determine if these coefficients 
are also universal, in the sense of (4). We have calcu¬ 
lated the corresponding ratios in the free scalar theory 
for n = 2, 3: (T 2 /C't = 2/9 and cts/Ct = 87r/(81-\/3) (ap¬ 
pendix B). It would be interesting to investigate these 
quantities in other theories, e.g., for a free fermion. 

While a{9) is only known implicitly for the free fields 
[5-7] and for holographic CFTs [12], the so-called ex¬ 
tensive mutual information model [36-38] produces the 
simple closed form expression 

aExt(6') = tt^Ct/S (1 -I- (tt- 6»)cot6») , (14) 

where we have fixed the overall normalization to recover 
(4). Interestingly, Ct arises from a pre-factor, which is 
the derivative of the scaling dimension of a twist operator 
with respect to the Renyi index. This derivative was 
recently shown to be tt^CtI^A: for general CFTs [46] - see 
appendix C for further comments. This simple expression 
(14) is also shown in Fig. 2 and exhibits the same nearly 
universal behavior as the other theories. 

Finally, we observe that the holographic result for 
a{9)/CT takes the smallest values in Fig. 2 across the 
whole range of 9. The fact that the results for all the 
other CFTs considered, both interacting and free, per¬ 
sistently lie above the holographic curve is certainly re¬ 
markable. It would be of interest to investigate whether 
the holographic result represents some kind of univer¬ 
sal lower bound. This is reminiscent of the conjecture 
that 77 /s is minimized by holographic CFTs dual to Ein¬ 
stein gravity [47]. In that case, however, the introduc¬ 
tion of higher curvature terms was found to invalidate 
the claim [48, 49], whereas ct/Ct remains unaffected by 
such terms. A challenge to this conjecture would arise 
from holographic models where the corresponding bulk 
surface is no longer the same as in Einstein gravity [21]. 
This would modify the functional form of a{9) and hence 
the lower bound might be violated for some values of the 
bulk couplings. 
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A more speculative suggestion comes from the obser¬ 
vation that holographic models describe CFTs at infinite 
coupling. It is tempting to guess that for general CFTs, 
the corner function should lie between those of the holo¬ 
graphic models and the free scalar field in Fig. 2. While 
the Wilson-Fisher CFTs are outliers at present, this pos¬ 
sibility can still be accommodated by the estimated error 
bars in Table I. Improved numerical results might soon 
clarify this situation [50]. 
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Appendix A: Holographic calculations 


In holographic theories dual to Einstein gravity, the EE is computed using the Ryu-Takayanagi prescription [18, 19]. 
This states that the EE of a region V in the d-dimensional boundary of an asymptotically AdS^+i spacetime is 
proportional to the area ^( 7 ) of the codimension-2 bulk surface 7 which is homologous to V in the boundary (in 
particular, dy = dV) and extremizes the area functional. That is, we have 


S{V) = ext 

7~y 


4G 


(Al) 


We are interested in the vacuum state of various three-dimensional holographic theories, so our bulk geometry will 
be four-dimensional Euclidean anti-de Sitter space, which we write in Poincare coordinates as 

72 

ds^ = dt^, + dp^ -I- p^dfj )^), (A2) 

where is the Euclidean time. This metric is a solution of Einstein gravity coupled to a negative cosmological 
constant. 


lo 


1 

IQttG 




(A3) 


provided we identify both length scales, L = L. We can compute the HEE for a region V delimited by two straight 
lines which intersect at the origin forming a corner with opening angle 9: V = {t^ = 0, p > 0, |(()| < 9/2}. The bulk 
surfaces 7 can be parameterized as 2: = ph{(j)), where h{(j)) approaches 5/p as <() —)■ ±d/2, with 5 being the UV cut-off. 
By symmetry, we also have ft.((/)) = h{—(j)) and hence d^h\^^Q = 0. The final result for the HEE, obtained using (Al), 
reads 


A 


i 

^ 5 


a^{9) ln(f/d)+0(l) 


(A4) 
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where we have introduced an IR regulator scale, = £, to ensure that the EE does not diverge. The function aE(d) 
is implicitly given by the following expressions [ 12 ] 


asiK) 

0 {ho) 


2G 


pOO 

/ dy 

1- J 

Jo 

L V 


dh 


+ /ig(l + y'^) 

2 + /i2(1 + j/2) 

2 hVI±E 


Vl + h?\/{hl - h‘^){hl + (1 + hDK^) 


(AS) 

(A 6 ) 


When bulk gravity theory contains higher curvature terms, the Ryu-Takayanagi prescription for HEE must be 
revised. In particular, the area functional in (Al) must be replaced by a new gravitational entropy functional (see 
e.g., [27-29]), just like the Bekenstein-Hawking formula for the black-hole entropy is replaced by Wald entropy [51-53] 
in that case. Schematically, we have 


S{V) = ext S'g„,( 7 ), 

■^r^V 


(A7) 


where the new functional depends on the details of the higher curvature action. The form of 5'grav for the action 
considered in the main text 


is given by [21, 32] 


I = 


= 


IGttG 


fd'^yVh 
. 4G 


“Pj -l- 7? -f Ij^X\RX^ X2X^ 


1 -f AiL'‘(A 4 -b 2RTI) + 4A2L®A’47^ 


(AS) 


(A9) 


When evaluated on empty AdS 4 , we find R = —12/L^, A 4 = 24:/L'^. Further, it is not difficult to show that the 
combination \/hTZ is a total derivative. In particular, one finds 


VhTZ = 


2(^-{l + 2h'^)d^h^-d4,h^ + {h + h^)dlh^ ^ 
pK^ (l + K^ + d'P 


dsh 


p sJl + K^ + 


(AlO) 


Using these results, one finds that the extremal bulk surface is unchanged and that a{6) is only modified by an overall 
factor, 


i{9) = aaE^d) with a = 1-|-24Ai-b O(A^). 


(All) 


Let us close this section by mentioning that we have actually computed a{0) and for the following broader class 
of higher-curvature theories - see [ 21 ] for details. 


/ = 


IGttG 


J d‘^Xy/g 


6 


— -b i? -b (Aii?^ -b X2 R^vR^^ + AGBA 4 ) 


(A12) 


{X^ qR^ -b Xi xRX^ -b LR {XiftR'^ -b X2pFR'X/^ -b Ao, 2 A’ 4 ) 


Note that the couplings Ai^ and Ao ,2 above correspond to Ai and A 2 in (A 8 ), respectively. The final expression of the 
corner coefficient and the corresponding charges a and Ct take the form 

a{9)=aaE{0), a = aaE, and Gt = q;Gt,e) (A13) 


where to leading order in the dimensionless couplings, the overall coefficient is given by 

a = l- 24Ai - 6 A 2 -b 432A3 ,o + 24Ai,i - 6912A4,o - 576A2,i -b ^(A^). 


Hence, we have 


a{9) aE(9) 


Ct Ct,e 

for all the boundary CFTs which are dual to (A12). 


and 


a (Te 


(AM) 


(A15) 
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Appendix B: Field theory calculations of a 

The first fourteen coefficients in the Taylor expansion of a{0) around 0 = tt were computed numerically for the 
cases of a free massless scalar and a free massless Dirac fermion using quantum field theory techniques in [5-7] . The 
first nonvanishing coefficients correspond to Uscaiar and CTfermion, which can be obtained by evaluating the following 
complicated integrals 


where 


^+oo ^+oo 

^scalar = —27r / dm / db ^ H a{l — a) m sech^ (jrb), 

Ji/2 Jo 

m cosech^ {-Kb) , 


/l/2 

f*+oo 


f*+oo 


<^fermion — 47r 


dm 


db 


ll /2 


F 


H = 

h = 

c = 




167ra(a — 1) ’ 

2 (a(a — 1) + sin^(7ra) 

(cos(27ra) + cos (ttVI — dm^)) ’ 

2 ^““^ 7 ra(l — a) sec (| ( 2 a + \/l — r(| — a + — dm^) 


= - 


mr(2 — a)2 r(a — | + \'J^ — dm^) 

r(-a) [tt sinh (^) + i cosh (^) (a + f + 5 ) - (« " T + ^))] 


22 o+i^r(a + 1) r(—a — ^ + 5) r(—a + ^ + 5) (cos(27ra) + cosh(7r/r)) 
Xi = “Ai” with a replaced by (1 — a), 

T = \/h{a^ — a + {h + l)m‘^), 

F 2 ’ 

(Sttc^ (to^ + 1 ) AiT + Stt/t (to^ + l) X 2 T - ch) - IGTra^T (c^Ai + hAz) 

a (—Sttc^to^AiT — STrhm^X 2 T + ch) + Sira'^T (c^Ai + / 1 A 2 ) — ch(h + 1 )to^ 
8c h (a^ — a + 


F 


F2 = 

fi = 

a = 


{ 2 a — l)fj, ’ 

\/ — 1 , 

f 6 + I for the scalar, 


i b 


for the fermion, 


(Bl) 

(B2) 


(B3) 


and denotes the digamma function i.e., ^ Inr(z). Notice that Eqs. (Bl) and (B2) look very different 

and without further insight, we find no reason to believe that these integrals should produce either simple or similar 
results. 

It is possible to compute integrals (Bl) and (B2) numerically with arbitrary precision (although, of course, the 
computation time increases considerably as we increase the precision). Our results indicate that both Eqs. (Bl) and 
(B2) exactly produce the results predicted assuming that cr/Cy is given by the universal constant ct/Ct = 7r^/24, i.e., 


^scalar = 7^ = 0.00390625 , 
256 


CTfermion = = 0.0078125 . 

iZo 


We have verified this result to a precision of approximately one part in 10^^. In particular, we find 
^scalar = 0.00390625000000(5) , CTfermion = 0.00781250000000(7), 


(B4) 

(B5) 


where the numbers in brackets are out of the range of accuracy of our computation. Let us note that the fact that the 
previous numerical results seemed to satisfy CTfermion = 2CTscaiar was observed in [5-7] , but no explanation was given. 
According to our conjecture, the reason comes simply from the well-known result that = 2(7®'^*^*^'' [17]. 

In fact, Renyi entropies contain regulator independent a„( 6 *) functions analogous to a{0), which corresponds to 
n = 1. In the regime of a nearly smooth entangling surface, these Renyi corner functions also behave as 

a„( 6 > tt) = CT„ (tt - 6 »)^ , (B 6 ) 















which straightforwardly generalizes the definition of cr to n ^ 1. The general expression we obtain for (T„, corresponding 
to the nth Renyi entropy (with n > 1) is given by [7] 




n—1 


E 


4fc(n 


1 ) 



dmm^/rn^ 


l/477/c/„ , 


(B7) 


where i7fc/n takes the same form as H in (B3) but with a = k/n instead. For n = 2 and 3, we find 


0-2 





l/4i7i/2 = 0.002140857992548703571747488816869325810507(4), 


1/4 (iJi /3 + H 2 / 3 ) = 0.001701632393277135955281871908373362473448(7). (B 8 ) 


These numerical results seem to fit the analytic following expressions 

= 108^ ’ 


(B9) 


within the numerical accuracy range - these analytic expressions were originally inferred with 20 significant digits in 
(B 8 ) and then tested by extending the accuracy to 40 significant digits. With these expressions, the corresponding 
ratios are given by 


(72 2 (T 3 87r 

Ct ~ 9" Ct~ 81^3 ■ 


(BIO) 


As noted in the main text, computing these quantities, and more generally OnjCT^ for other theories and investigating 
whether the results in (BIO) are universal, would be of great interest. 


Appendix C: Corner coefficient a(ff) for the extensive mutual information model 


In this section, we outline the derivation of the corner coefficient aBxt(^)) Eq. (C 8 ), in the so-called extensive mutual 
information model of [36, 38]. As its name suggests, this model is characterized by the special property that the 
mutual information I{A,B) satisfies the extensivity property: I{A,B U C) = I{A,B) + I{A,C). We note that the 
expression for aExt(0) first appeared in Ref. [36], and was later independently obtained by Swingle [37]. Here, we shall 
follow Swingle’s heuristic presentation in terms of twist operators, but emphasize that the end result is the same. 
While this model can be extended to any number of spacetime dimensions, we only consider d = 3. A key role is 
played by the twist operator Kn [H], which is a line-operator that introduces a branch cut at the boundary dV in the 
n-fold replicated theory. The expectation value of Kn in the n-replicated theory yields {Kn\V]) = Tr(py), as required 
to calculate the nth Renyi entropy. Inspired by results for twist operators in d = 2 CFTs, Swingle made the following 
ansatz for the twist field in higher dimensions: 


Kn\y] = exp 



(Cl) 


where n is the (spatial) unit vector normal to the boundary, and ^ is a two-component vector field (in d = 3) defined 
on dV. For simplicity, (f is assumed to be a Gaussian field, and thus entirely determined by its two-point function. 
The latter can be taken to be 


((/*(f)</>^(0)) 


bi d'-J -I- 62 


(C2) 


where bi are real coefficients and x’’ denotes the components of the unit vectors x = x/\x\. In contrast to [37], we 
shall keep track of how the coefficients, 61 and 62 , enter in the final answer for the EE and we have also extended the 
correlator (C2) to include the 62 term. The EE is given by S'(F) = —^n{Kn\V])\n=l^ which simplifies to the following 
double integral over dV 


S{V) 


1 

2 


dn{al) 


n=l 


JdV JdV 


{(t>\Xi) fP {X2)) , 


(C3) 
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by virtue of (Cl) and the Gaussianity assumption. 

Specializing to the case of a wedge shaped region with opening angle 0, the integral can be explicitly evaluated and 
yields a result of the expected form 


S{V) =B-- a( 6 ») In {e/6) + 0{S/e)° , (C4) 

where <5 is a UV cutoff, i an IR scale, and where the corner coefficient reads 

a{d) = ^ dnial)\^^^ (1 + (tt - 0) cot 9 ), (C5) 

which holds for 0 < 0 < 2tt. Hence both the bi and 62 terms in the correlator (C2) make the same contribution to 
the logarithmic term up to an overall factor. 

Following [37], we interpret as the scaling dimension of the twist operator and hence we can apply the recent 
result of [46] 


9n{an)\n=l ^ ’ 

which holds for general CFTs in d = 3 - see also [54]. Substituting this expression into (C5), we set 

26 i +62 = 6 /tt 

to produce the correct normalization for 9 near tt, i.e., a{9) ~ tt^Ct (tt — 0)^/24. Hence we find 


(C 6 ) 


(C7) 


aExt(^') = (1 + (tt - 0 ) cot d) , 


(C 8 ) 


which is plotted in Fig.2 of the main text. We also note that, as previously observed in [36, 37], the above corner 
coefficient diverges as ^ k/ 0 for 0 —)■ O"*". The coefficient can be easily determined: 


K = 



(C9) 


which yields a distinct ratio for k/Ct compared to the free scalar and Dirac fermion field theories, as well as to the 
holographic CFTs. Further, it can be checked that a{9) satisfies the non-trivial inequalities: 

0 ( 0 ) >0, a'(0)<O, a'(0)+sin(6»)a"(0) > 0, (CIO) 

obtained using strong subadditivity of entanglement and Lorentz invariance [5]. 

It is not clear at present why the simple extensive mutual information model, when interpreted in terms of twist 
operators, captures the correct factor of in the d —>■ tt limit. One non-universal property of this model is that 
will appear as an overall pre-factor of the EE associated with any region V. This includes the case where H is a disk, 
which in turn implies that Ct fixes the RG monotone F. In this case, the EE is easily obtained as: 

S{disk) = B^-F^^, + 0{d/R) (CIl) 

d 

F^^, = ^{2h+b2)dn{al)l^^, (CI2) 

where R is the radius of the disk. Ref. [37] found the same expansion for i? ^ d, but did not specify the form of 
the regulator independent constant F. Here we not only identify its sign, but also its value in terms of the model 
parameters. In particular, we see that Eexi contains the same combination (26i -|- 62 ) that appeared above in aExt(0)- 
Hence using the same values for that linear combination and the n-derivative, we obtain: 

Eex* = ^ , (CI3) 

o 


which is consistent with the fact that F > 0 for a GET. We note that for a CFT holographically dual to pure Einstein 
gravity, F^ = Of course, the latter differs from (C13) but as illustrated in the main text, there is no universal 

relation between F and Ct- Nevertheless, CETs described by the extensive mutual information model [59] and those 
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dual to Einstein gravity seem to share the non-generic feature that Ct controls the EE for any region and so it would 
be interesting to see whether a deeper connection exists between this model and holography. 
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